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ABSTRACT 


DEFINITION OF SYMBOLS 


In free piston Stirling engine research the 
integrity of both amplitude and phase of the dynamic 
pressure measurements is critical to the 
characterization of cycle dynamics and 
thermodynamics. It is therefore necessary to 
appreciate all possible sources of signal distortion 
when designing pressure measurement systems for this 
type of research. This paper discusses the signal 
distortion inherent to pressure transmission lines. 
Based on results from classical analysis, guidelines 
are formulated to describe the dynamic response 
properties of a volume -terminated transmission tube 
for applications involving helium-charged free 
piston Stirling engines. The scope and limitations 
of the dynamic response analysis are considered. 


INTRODUCTION 

It is often impossible to flush-mount a dynamic 
pressure transducer in the expansion space or the 
compression space of a free-piston Stirling engine. 
In some cases, direct exposure to the severe gas 
temperatures would threaten to destroy the 
transducer; In other cases, the engine design 
precludes straightforward transducer access. When 
flush-mounting is not possible, the transducer is 
physically removed from the working space via a 
volume-terminated transmission tube (some 
termination volume at the instrument end of the tube 
is an inevitable consequence of transducer 
installation. ) 


The following symbols are defined in order to 
facilitate the discussion of pressure transmission 
lines and their use in Stirling engine research: 


SYMBOL 

DEFINITION 


FPSE 

free piston Stirling engine 

L 

tube length 

[ cm ] 

D 

tube inner diameter 

t cm ] 

V t 

tube volume 

[ cm 3 ] 

Vi 

instrument volume 

( cm 3 ] 

Vr 

volume ratio « V t / V t 


T * 

gas mean temperature 

[ K ] 

P * 

gas mean pressure 

[ bar ] 


fluid viscosity 

[g/cm* s 

V 

kinematic viscosity 

[cm 2 / s ] 


specific gas constant 


7 

ratio of specific heats 


c 

gas speed of sound 

- V R, 7 t 9 


w 

angular frequency 

[ l/s ] 

f 0 

resonant frequency of 



tube and volume system 

t Hz ] 


In these situations, the instrumentation engineer 
must ensure that errors introduced in the pressure 
signal by the transmission tube remain within 
acceptable tolerances. In Stirling engine research, 
the error tolerances on dynamic pressure 
measurements are particularly stringent since small 
errors in either amplitude or phase of the 
measurement can cause significant errors in 
important calculations, such as the PV power. 

This paper will develop some useful guidelines for 
the design of pressure transmission lines for use in 
helium-charged free piston engine applications. 


f ang Stirling engine's operating 

frequency ( Hz ] 

<p pressure phase shift [ deg ] 

z Iberall ' s damping parameter 

X Iberall 's attenuation parameter 

B&T Bergh and Tij deman response 

analysis 


DYNAMIC RESPONSE OF A VOLUME -TERMINATED TUBE 

The general analysis of oscillatory pressure 
transmission through a volume-terminated tube is a 
classical problem in instrumentation engineering. 
Since a transmission tube will distort both the 
amplitude and phase of a dynamic pressure 
measurement, it is necessary to perform a complete 
dynamic-response analysis to theoretically bound the 
measurement errors. 

To date, the most comprehensive analytical treatment 
of this problem has been presented by Iberall[l]. 
In his 1950 treatise, Iberall solved the fundamental 
equations of fluid flow in order to establish how an 
oscillatory pressure is transmitted through a tube 
and volume system. The resulting dynamic response 
curves depict both the amplitude distortion and the 
phase shift as a function of oscillation frequency. 


There are a number of available computer algorithms 
based on Iberall* s results which have been 
successfully validated by laboratory experiments. 
One example is the code presented by Goldschmied [2] . 
Another is the algorithm developed at the NASA Lewis 
Research Center by Nyland, Englund and Anderson[3]; 
the NASA code implements a method presented by Bergh 
and Tijdeman[4] which extends Iberall 's mathematics 
to handle more complicated tube and volume 
geometries . 

In the Stirling Engine Branch at NASA Lewis, the 
Bergh and Tij deman code has become a standard tool 
in the evaluation of pressure instrumentation 
systems. This paper will draw upon computer results 
from the B&T code along with analytical results from 
Iberall’s paper to demonstrate basic properties of 
pressure transmission response. 


Figure 1 represents a typical two-part response 
curve derived from Iberall’s mathematical technique: 



Figure Is TYPICAL DYNAMIC RESPONSE BEHAVIOR FOR PRESSURE 
TRANSMISSION THROUGH A VOLUME-TERMINATED TUBE 


It can be seen that for frequencies which are low 
with respect to the primary resonant frequency, the 
dynamic response curves are reasonably flat; dynamic 
pressures at these frequencies will be transmitted 
predictably with little distortion. It is necessary 
to take advantage of this flat portion of the 
response curve when designing a pressure 
transmission system. The instrumentation engineer 
needs to be equipped with basic analytical tools in 
order to perform such a design. 


DESCRIPTION OF THE UNDERDAMPED FLOW REGIME 


In general, the response curves for a tube and 
volume system will be influenced by tube length, 
tube diameter, volume ratio and gas properties 
(which are a function of gas mean temperature and 
gas mean pressure.) However, an evaluation of 
pressure transmission lines in existing helium- 
charged free-piston research engines has yielded 
some interesting observations: although the 

resonant frequencies of these lines are strongly 
affected by tube length, volume ratio and gas mean 
temperature, they are essentially independent of 
tube diameter and gas mean pressure. 


Using results from the Iberall analysis, it will be 
shown that this behavior is characteristic of 
transmission tube response when the working fluid in 
the tube is ’underdamped*. First, the conditions 
for underdamped flow will be discussed; and it will 
be shown that the helium in pressure transmission 
lines of typical FPSE applications will behave as an 
underdamped fluid. 

The dimensionless parameter used by Iberall [1] to 
indicate the amount of damping in a fluid flow is z : 

D 2 W D 2 IT f 

z = = [ eq. 1 ] 

4 V 2 1/ 

Here, V is the kinematic gas viscosity, and f is the 
frequency of the pressure oscillation. The fluid is 
considered underdamped when z > 100. 


For a given FPSE application, f. is fixed by the 
engine’s operating frequency, and V is a fixed 
function of gas mean temperature and gas mean 
pressure; it then follows that the tube diameter, D, 
is the only independent variable in this equation. 
If z is replaced with the threshold value of 100, 
equation 1 can be rearranged to give an expression 
for the minimum tube I.D. which will guarantee 
underdamped gas flow in the tube; 


D 


min ~ 


200 


eq. 2 ] 


2 




Table I calculates D ain for a helium-charged engine 
with f -ng « 100Hz at various values of gas mean 
temperature and pressure. 


x g ‘g " »in 

[ K ] [ bar ) [ cm 2 /s ] [ cm ] 


300 

150 

.0088 

.075 

750 

150 

.0390 

.157 

1050 

150 

.0676 

.207 

300 

100 

.0129 

.091 

750 

100 

.0585 

.193 

300 

50 

.0628 

.200 


TABLE I: D Bin §100Hz vs. T g and P g for Helium 
( v values from gas tables ) 


It can be seen that for many conventional tube sizes 
(.32cm, .64cm etc.), D tuba will exceed D Bin in 100Hz 
FPSE applications; thus the helium is underdamped. 

This paper will argue that underdamped flow is not 
only achievable in most FPSE applications, it is 
also desirable. The analysis of transmission tubes 
becomes simplified for low damping; and more 
important, the pressure phase shifts are 
considerably less severe when the gas flow is 
underdamped . 

In consequence, the instrumentation engineer should 
select tubes for a pressure transmission line which 
ensure underdamped helium flow at the specific 
application frequency. Equation 2 provides the 
necessary guidance. 


THE PRIMARY RESONANCE FREQUENCY 

As previously illustrated in figure 1, the most 
important design consideration for a pressure 
transmission system is the primary resonant 
frequency. For a complete description of resonance 
phenomena in a volume-terminated transmission tube, 
the reader is referred to Iberall’s complete 
dissertation[l] ; The basic results developed for 
the underdamped flow regime, however, represent a 
significant simplification of the generalized 
response analysis. Two of these simplified results 
can be summarized as follows: 

1) The primary resonant frequency for a transmission 
tube at a fixed volume ratio will vary in proportion 
to the quantity C/L. Here, C represents the gas 
speed of sound and L is the tube length. 

2) For a fixed value of C/L, the frequency of 
primary resonance is solely a function of the volume 
ratio. The theory does not provide a simple 


analytical expression for this function; for 
convenience, it will be denoted by G(Vr) in this 
paper. 

A combination of these results implies that the 
primary resonant frequency can be described by the 
following proportionality: 

C 

f 0 « G( Vr ) ( eq. 3 ] 

L 

Or, since the gas speed of sound, C, is proportional 
to the square root of the gas mean temperature: 



L 


This represents a generalized expression for the 
primary response resonance when the fluid flow in 
the transmission tube is underdamped. In this 
formula there are no explicit dependencies on tube 
diameter or gas mean pressure; a result which is 
consistent with the observations mentioned earlier 
for transmission lines used in helium-charged 
FPSE’s. 

A specific equation for underdamped helium systems 
can be developed in the form of equation 4. 
Resonant frequency information derived using the 
Bergh and Tij deman code for a family of such systems 
is presented in Table II ( the tube length of the 
systems is 1 cm, and the helium temperature is 
300K. ) Primary resonance is calculated as a 
function of volume ratio. Using this information, 
the helium-specific equation will be formulated 
empirically. 


VOLUME 

RATIO 

fo 

[ Hz ] 

.01 

2.583E4 

.1 

2.258E4 

.5 

1.537E4 

1.0 

1.181E4 

2.0 

8.739E3 

3 . 0 

7. 247E3 

5.0 

5.687E3 

7.5 

4.673E3 

10.0 

4.058E3 

12.5 

3.713E3 

15.0 

3.326E3 


TABLE II: B&T PREDICTIONS FOR f 0 IN AN UNDERDAMPED He 

TRANSMISSION TUBE (T -300K; P -150bar; L-lcm; 
D- . 325cm; 2-1878 @ 100Hz) 


Figure 2 is a plot of the computer-predicted 
resonant frequencies from Table II as a function of 
volume ratio. Superimposed on the plot is a curve 
that was empirically derived to fit the points. 
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20784.8830 ( Vr* 6574 + .7365)’ 1 


PRESSURE AMPLITUDE DISTORTION 



0 5 10 15 

VOLUME RATIO 

Figure 2: f 0 vs Vr FOR A 1 cm TUBE WITH UNDERDAMPED HELIUM 
^ T g «300K 

An expression for the primary resonant frequency of 
any underdamped helium transmission tube can now be 
deduced by forcing the curve-fit equation into the 
mathematical form of equation 4: 

f 0 (Hz) ^ 20784.8830 ( Vr- 6574 + . 7365 )* x 

= 1200 v' 300 K 1 


1 cm ( Vr* 6574 + . 7365 ) 


- 1200 n/ T g (K) 1 1 

L (cm) ( Vr* 6574 + . 7365 ) 

[ eq. 5 ] 

Table III lists the resonant frequencies calculated 
using this equation alongside the B&T predictions 
for several underdamped helium systems. Good 
correlation exists between predictions of equation 
5 and the computer results. 


L 

C cm ] 

D 

( cm ) 

Vr 

T* 

( K ] 

P. 

[ bar ] 

1, : Hz 
( B&T ] 

t, : Hz 
[ EO.5 ] 

DEVIATION 

1.0 

.16 

.25 

300 

100 

18650 

18256 

-2. IX 

1.5 

.16 

2.5 

300 

150 

5250 

5406 

3.0% 

2.0 

.32 

0.5 

300 

200 

7748 

7583 

-2.1% 

2.5 

.32 

5.0 

500 

150 

2888 

2967 

2.7% 

3.0 

.64 

.75 

750 

150 

6830 

7003 

2.5% 

5.0 

.64 

7.5 

1050 

150 

1701 

1729 

1.6% 

10.0 

.64 

10.0 

300 

150 

405 

394 

-2.7% 


TABLE III: COMPARISON OF EQUATION 5 WITH B&T PREDICTIONS FOR 
7 TRANSMISSION TUBES WITH UNDER0AMPED HELIUM 


Equation 5 is therefore a legitimate predictor of f 0 
in pressure transmission tubes which incorporate 
underdamped helium. Since it explicitly includes a 
volume-ratio dependence, it will be more generally 
valid than the conventional * quarter wave ’ equation 
for acoustic resonance of a closed-ended tube. 


In engineering problems which can be analyzed as 
second order systems, a knowledge of the first 
resonance will be sufficient to characterize 
amplitude response at frequencies that are low with 
respect to f 0 . In other words, although f 0 will 
shift as system parameters are varied, the amplitude 
distortion at some low frequency, f, will remain a 
fixed function of f/f 0 . 

Frequency response information for a family of 
underdamped helium systems is presented in Table IV. 
Entries in the table indicate that the tube and 
volume systems under consideration exhibit this 
second-order property of amplitude response. 


Vr 

fo 

[ Hz ] 

f 10% 

C Hz ] 

f 10% 
fo 

f5X 
C Hz ] 

f5X 

fix 

[ Hz ] 

fix 

fo 

.01 

2.583E4 

7.056E3 

.273 

5.085E3 

.197 

2.308E3 

.089 

.1 

2.253E4 

6.239E3 

.276 

4.493E3 

.199 

2.040E3 

.090 

.5 

1.537E4 

4.396E3 

.286 

3.169E3 

.206 

1.440E3 

.094 

1.0 

1.181E4 

3.450E3 

.292 

2.485E3 

.210 

1.133E3 

.096 

2.0 

8.739E3 

2.588E3 

.296 

1.870E3 

.214 

8.404E2 

.097 

3.0 

7.247E3 

2.157E3 

.298 

1.559E3 

.215 

7.103E2 

.098 

5.0 

5.687E3 

1.701E3 

.299 

1.231E3 

.216 

5.586E2 

.098 

7.5 

4.673E3 

1.404E3 

.300 

1.012E3 

.217 

4.599E2 

.098 

10.0 

4.058E3 

1.221E3 

.301 

8.802E2 

.217 

4.000E2 

.099 

12.5 

3.713E3 

1.118E3 

.301 

8.046E2 

.217 

3.673E2 

.099 

15.0 

3.326E3 

9.977E2 

.300 

7.23 1E2 

.217 

3.287E2 

.099 


TABLE IV: B&T FREQUENCY RESPONSE PREDICTIONS FOR A PRESSURE 
TRANSMISSION TUBE WITH UNDER0AMPED HELIUM ( T - 
300K; P -ISObar; L-lcm; D-. 325cm; 2-18780100H! ) 


The columns labelled f 10% , f 5% and f 1% in the table 
represent the frequencies at which the amplitude 
ratio becomes 1.10, 1.05 and 1.01 respectively. 
Independent of the volume ratio, these frequencies 
are essentially fixed percentages of the resonant 
frequency. Specifically: 

f i% 2 .10 f 0 [ eq. 6 ] 

f 5% ? .21 f 0 [ eq. 7 ] 

f 10% ? .30 f 0 ( eq. 8 ] 

Dynamic pressure signals oscillating at frequencies 
below f 1% will pass the transmission line with an 
amplitude distortion of less than 1Z. Similar 
definitions hold for f 5% and f l0% . Moreover, linear 
interpolation can be used to approximate 
intermediate values of amplitude distortion. For 
example, f 3% can be interpolated between f 1% and f 5% : 

f 3% - + %(f 5 % “ *i%) - -15 f 0 [ eq. 9 ) 

Thus, after calculating the resonant frequency for 
a given pressure measurement system, error bands on 
the amplitude distortion are easily established. 

Consider a hypothetical problem from the 
instrumentation designer's point of view: It is 
necessary to measure the amplitude of dynamic 
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pressures in a 100 Hz engine to within 3 Z ; Then, f 3% 
is taken to be 100 Hz. Based on the result of 
equation 9, the transmission lines must be designed 
to resonate at f 0 > f 3% * .15 * 667 Hz. The higher 
f 0 can be made, the better in terms of pressure 
amplitude integrity. 


PRESSURE PHASE DISTORTION 

In second order systems that are completely 
undamped, there will be no phase shift at 
frequencies below the primary resonance. In 
underdamped systems, the phase shifts can generally 
be made small by careful system design. This result 
is particularly significant to FPSE research, where 
it is often necessary to resolve the pressure phases 
to within one tenth of one degree. 

Unfortunately, the detailed analysis of phase 
distortion is somewhat more complicated than that of 
amplitude distortion. The phase shift will depend 
not only the factor f/f 0 , but also on the amount of 
damping present in the transmission system. ( In 
the case of amplitude distortion, as long as the 
system is underdamped, the actual amount of damping 
is unimportant to low-frequency response. ) 

The appendix offers a heuristic development of phase 
distortion behavior in underdamped systems. The 
resulting approximation describes the functionality 
of ^ for a fixed volume-ratio: 

L 2 f 1 * 5 

^ a [ eq. 10 ] 

D /"p7 

Table V presents Bergh and Tijdeman predictions for 
the phase distortion in a specific underdamped 
helium transmission tube. ( L is 1 cm, D is .325 cm 
and P g is 150 bar.) The phase shifts for a pressure 
oscillating at 100 Hz are tabulated as a function of 
volume ratio. 


VOLUME RATIO <p § 100 Hz 


.01 

-.00006 

0.1 

-.00007 

0.5 

-.00012 

1.0 

-.00019 

2.0 

- . 00027' 

3.0 

-.00042' 

5.0 

-.00064' 

7.5 

-.00093' 

10.0 

-.00121' 

12.5 

-.00144' 

15.0 

-.00179' 


TABLE V: B&T PREDICTIONS FOR <t> 0100Hz IN AN UNOERDAMPED 

HELIUM TRANSMISSION TUBE (P -150bar; T -300K; 
L-lcro; D-. 325cm; 2-1878 ) 9 9 


Using the scaling rules of equation 10, information 
in Table V can be applied generally to any 
underdamped helium system as follows: 

(L/lcm) 2 (f/lOOHz) 1 ' 5 

^general — Arabia 

(D/ . 325cm) V P g /150bar 
L 2 f 1 ’ 5 

3.98E-3 — - * tabl . 

D v' P g 

[ eq. 11 ] 

Although the absolute values of <p in Table V are 
orders of magnitude less than .1°, it can be seen 
that an imprudent system design will result in 
unacceptably large phase shifts. Consider, for 
example, a 150 bar engine operating at 100 Hz. A 
transmission tube with L=10cm, D=.16cm and Vr=10 can 
distort the phase of the pressure signal by as much 
as .25°; this by virtue of equation 11 and Table V. 

The key point here is that phase shifts in 
underdamped systems are generally small; however, 
increasing the tube length of the system will 
magnify the phase distortion by a factor of L 2 . 


LIMITS OF APPLICABILITY 

The specific guidelines derived in this paper apply 
exclusively to volume- terminated transmission tubes 
which incorporate underdamped helium as the internal 
fluid. There are, moreover, underlying assumptions 
built into Iberall’s fundamental mathematics which 
further limit the applicability of the results. 

As one notable example, the fundamental flow 
equations at the root of the analysis have been 
linearized based on an assumption that dynamic 
pressure amplitudes are small compared to the mean 
pressure. Fortunately, Iberall addresses the 
implications of this assumption for situations where 
the pressure amplitude is large with respect to mean 
pressure. Specifically, the fundamental of a large 
pressure oscillation will be transmitted in 
accordance with the 1 inear ized-mode 1 predictions; 
however, higher-order harmonic distortion may be 
artificially induced in the pressure measurement. 
Therefore, in FPSE research, where dynamic pressure 
amplitudes as large as 202 of the mean pressure are 
not uncommon, harmonic distortion in measured 
pressure signals may be a phantom of the 
transmission line. 

Another shortcoming of the theory is its inability 
to account for complicated thermal interactions 
between the gas and the tube. If, for example, 
dynamic pressure in a FPSE’s high-temperature 
expansion space is being measured by a water-cooled 
transducer at the end of a transmission tube, then 
severe temperature gradients will be induced along 
the tube. The extent of resulting anomalies in the 
pressure measurement cannot be predicted using 
methods based on Iberall’s classical analysis. 
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CONCLUDING REMARKS 

At the NASA Lewis Research Center, there is an 
effort in progress to examine a more generalized 
solution to the tube and volume response problem 
using state of the art computational fluid dynamics 
methods. In conjunction, an instrumentation 
research program is under development to provide 
experimental verification of the improved analysis. 

Until results from this effort are available, the 
classical analysis should be applied with 
conservatism. In practice, the Iberall theory can 
provide estimated error bands for amplitude 
distortion and phase shift; however, it is not 
adequate for calculating reliable correction factors 
to compensate for errors in the pressure 
measurement - 

There is one conclusion that has been implicit 
throughout this paper: Flush-mounting a pressure 
transducer, whenever possible, is preferred. This 
circumvents transmission line distortion, and it 
adds no unwanted dead-volume into the FPSE’s working 
space . 

When a flush-mount is not possible, then careful 
design of a pressure transmission system is 
imperative. This paper has developed some basic 
guidelines for instrumentation engineers to follow 
in the early design stages of a pressure measurement 
system. To summarize in general: the tube diameter 
should be large enough to ensure underdamped helium 
flow; and both volume-ratio and tube length should 
be minimized in order to maximize f 0 with respect to 
^•ng * 

The guidelines are approximate, and should be viewed 
as a starting point; they are no substitute for a 
full dynamic- response analysis. The computer codes 
cited in this paper are highly recommended; If 
properly used, they should prove invaluable to the 
instrumentation designer. Hopefully, the basic 
properties of dynamic response described in this 
paper will facilitate an efficient use of the 
analytical codes, and will aid in the effective 
design of pressure measurement systems. In free 
piston Stirling engine research, the importance of 
maintaining rigorous standards in the development of 
instrumentation systems cannot be overstated. 


APPENDIX: PHASE DISTORTION FOR UNDERDAMPED FLOW 

Two factors contribute to phase distortion in an 
underdamped transmission system at low frequencies 
with respect to f 0 . This appendix will examine these 
factors for systems with fixed volume-rat ios . 

One factor will vary as a function of f/f 0 ; an effect 
characteristic to both phase and amplitude 
distortion. For phase, the function is nearly 
linear in the low frequency region. In this 
appendix, the first factor will be denoted by X. 
Using equation 4 with a fixed volume-ratio: 


f f L 

X a a 

*0 /v 

The second factor represents the effect of damping 
on the phase shift. This factor will vary as the 
square root of the Iberall attenuation parameter 

jrm« 

sTT a J Ji ~ 9 • ( L / D ) 

( Through mathematical manipulation, the parameter 
X can be expressed in terms of 1 (z ; thus, as 
expected, high values of z express conditions of 
low damping. This is left as an exercise for the 
interested reader. ) 

Next, H will be replaced by the approximation /?*T g n 
( where p and n are constants, and n is close to, 
but less than 1. ) Then, combining the two phase- 
distortion factors leads to: 

<j> a X • v / - x 



L 2 f 1 * 5 


d 

Although this result is approximate, it provides an 
extremely useful description of the phase 
sensitivity to transmission tube parameters. 
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